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Abstract 

The free energy of a nonabelian gauge theory at high temperature T can be calcu- 
lated to order using resummed perturbation theory, but the method breaks down 
at order . A new method is developed for calculating the free energy to arbitrarily 
high accuracy in the high temperature limit. It involves the construction of a sequence 
of two effective field theories by first integrating out the momentum scale T and then 
integrating out the momentum scale gT. The free energy decomposes into the sum of 
three terms, corresponding to the momentum scales T, gT, and g'^T. The first term 
can be calculated as a perturbation series in g^{T), where g{T) is the running coupling 
constant. The second term in the free energy can be calculated as a perturbation series 
in g{T), beginning at order g'^. The third term can also be expressed as a series in g{T) 
beginning at order g^, with coefficients that can be calculated using lattice simulations 
of 3-dimensional QCD. Leading logarithms of T /{gT) and of gT /{g'^T) can be summed 
up using renormalization group equations. 



The perturbative infrared catastrophe of a nonabehan gauge theory at high temperature 
is one of the most important unsolved problems in thermal field theory. The problem, which 
was first identified by Linde in 1979 0, is that perturbative calculations of the free energy 
of a nonabelian gauge theory at high temperature T with weak coupling g break down at 
order g^. In contrast, there seems to be no obstacle to calculating the free energy of an abelian 
gauge theory like QED to arbitrarily high order in the coupling constant e. This problem 
casts a shadow over all applications of QCD perturbation theory to describe a quark-gluon 
plasma at high temperature, because the catastrophe afflicts any observable at sufficiently 
high order in g. The catastrophe also arises in applications of the electroweak gauge theory 
and grand-unified theories to describe phase transitions in the early universe, since they also 
involve nonabelian gauge theories at high temperature. The catastrophe can be avoided for 
static observables in hot QCD by calculating them directly using lattice simulations, but 
the computational resources that are required increase rapidly with T. Moreover, there are 
many problems, such as the calculation of dynamical observables or the calculation of static 
observables at nonzero baryon density, which are more easily addressed by perturbative 
methods than by lattice simulations. For these reasons, it is important to solve the problem 
of the perturbative infrared catastrophe. 

In this Letter, we solve the problem by constructing a sequence of two effective field 
theories. The first effective theory, which is obtained by integrating out the momentum 
scale T, is equivalent to thermal QCD at length scales of order l/{gT) or larger. The 
second effective theory, which is obtained by integrating out the scale gT from the first 
effective theory, is equivalent to thermal QCD at length scales of order l/{g'^T) or larger. 
If T is sufficiently large, the parameters of the two effective theories can be calculated as 
perturbation series in the running coupling constant g{T). The second effective theory is 
inherently nonperturbative, so that the effects from the scale g^T must be calculated by 
lattice simulations. 

To explain the origin of the infrared catastrophe, we first consider the case of the abelian 
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gauge theory QED. The leading term in the free energy is that of an ideal gas of electrons, 
positrons, and photons. The first correction due to the electromagnetic interaction is of 
order e^, and can be calculated using ordinary perturbation theory. Beyond that order, 
infrared divergences arise in perturbation theory and it is necessary to sum infinite classes 
of diagrams. The order correction to the free energy can be calculated by resumming 
"ring diagrams" , as in the classic calculation by Gell-Mann and Brueckner for a degenerate 
electron gas 0. This method was used to calculate the order correction to the free 
energy of a degenerate electron gas in the relativistic limit 0]. An equivalent resummation 
method was used to calculate the free energy to order for a massless scalar field with 
a 0^ interaction in the high temperature limit 0. More recently, this method was used to 
calculate the corrections to the free energy for hot QED of orders and [Q. The only 
apparent obstacle to calculating to still higher order in e is the increasing complexity of the 
sums and integrals that are encountered. 

We now consider the case of hot QCD. The resummation of ring diagrams allows the 
free energy to be calculated to order g^. The free energy for pure-glue QCD was recently 
calculated to order g^ by Arnold and Zhai using an equivalent resummation method [rU|. 
However, at order g^, there are contributions to the free energy that can only be calcu- 
lated using nonperturbative methods. To see where they arise, recall that static quantities 
in a thermal field theory can be calculated using the imaginary-time formalism, in which 
the propagator for a massless field is l/{uj^ + k^), where Un takes on the discrete values 
LOn = 2mTT for bosons and Un = (2?7, + 1)ttT for fermions. The only modes that can prop- 
agate over distances much larger than 1/T are the n = modes of the bosons. These 
are therefore the only modes that can give rise to infrared divergences that cause a break- 
down of perturbation theory. With the resummation of ring diagrams, the chromoelectric 
zero modes ( "electrostatic gluons" ) acquire a mass of order gT, while chromomagnetic zero 
modes ( "magnetostatic gluons") remain massless. When restricted to the chromomagnetic 
zero modes, the action of QCD reduces to that of 3-dimensional Euclidean QCD with cou- 
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pling constant g^T. The perturbation expansion for this theory is hopelessly afflicted with 
infrared divergences, but it is well-behaved nonperturbatively with a mass gap of order g'^T. 
Using simple dimensional analysis, one can show that the contribution to the free energy 
density from these modes is {g'^T)^T multiplied by a coefflcient that can only be obtained 
by a nonperturbative calculation. There is no analogous contribution in QED, because the 
restriction to magnetostatic zero modes gives a free field theory. 

We present the solution to the perturbative infrared catastrophe in the context of ther- 
mal QCD, although the solution applies equally well to any nonabelian gauge theory. Our 
starting point is the imaginary-time formalism for thermal QCD, in which static quantities 
are calculated in 4-dimensional Euclidean QCD, with the periodic Euclidean time r having 
period /3 = l/T. The partition function is 



^QCD 



(T) 



/ VA^^ exp J^^ dr J d'x^G%G^^}j , (1) 



where G^^, is the nonabelian field strength constructed out of the gauge field with coupling 
constant g. The free energy density is F = —TlogZ/V, where V is the volume of space. It 
can be calculated at any temperature T using lattice simulations of the functional integral. 
If T is large enough that the running coupling constant g{T) is small, one might expect 
to be able to calculate F{T) as a power series in g'^{T) using ordinary perturbation theory. 
However this perturbation expansion breaks down due to infrared divergences associated with 
the exchange of static gluons. The divergences due to electrostatic gluons can be eliminated 
by the resummation of ring diagrams, but a more elegant solution is to construct an effective 
field theory which reproduces the static correlation functions of thermal QCD at distances of 
order l/(gT) or larger. This theory, which we call electrostatic QCD (EQCD), contains an 
electrostatic gauge field Aq{-x) and a magnetostatic gauge field A^(x). Up to normalizations, 
they can be identified with the zero- frequency modes of the gluon field ^^(x, r) for thermal 



QCD in a static gauge |TT[. The lagrangian for this effective field theory is 

^EQCD = \Gt,Gt, + ^{D,AonD,A,r + ^miASA^ + 5£eqcd, (2) 



where Gij is the magnetostatic field strength with couphng constant g^. The lagrangian (0) 
has an SU{3) gauge symmetry. If the fields Aq and Ai are assigned the scaling dimension 1/2, 
then the operators shown explicitly in (^) have dimensions 3, 3, and 1. The term SCeqcd in 
(0) contains all possible local gauge-invariant operators of dimension 2 and higher that can 
be constructed out of Aq and A^. 

Renormalization theory tells us that static correlation functions at distances R ^ 1/T 
in the full theory can be reproduced in the effective theory with any accuracy desired by 
tuning the coupling constant gs, the mass parameter m^i, and the parameters in 5£eqcd as 



functions of T and the ultraviolet cutoff A of the effective theory [12]. The A-dependence of 
the parameters is cancelled by the A-dependence of the loop integrals in the effective theory. 
The coefficients of the operators in the effective lagrangian (0) for EQCD can be computed as 
perturbation series in g'^{T). At leading order in g, we have (73 = g\/T and m^i = Ncg'^T'^/3. 
The coefficients of some of the higher- dimension operators were recently computed to leading 
order by Chapman [jl^. Beyond leading order, it is particularly convenient to use a scale- 
invariant regularization scheme, such as dimensional regularization, which automatically 
eliminates all power ultraviolet divergences. In this case, the coefficient of an operator 
of dimension d is T^"*^ multiplied by a power series in g'^{T), with coefficients that are 
polynomials in log(T/A). The parameters satisfy simple renormalization group equations 
that can be used to sum up leading logarithms of T/A. 

The effective field theory EQCD can be used to calculate the partition function of thermal 
QCD, as well as static correlation functions at distance scales of order 1/ (gT) or larger. The 
partition function is 

^qcd(T) = e-f^^-^^^^'"" f'^^VA^.VAIexp [- j rf^x^EQCo) • (3) 

In the exponential prefactor, JqcbT^ can be interpreted as the coefficient of the unit operator 
which was omitted from the effective lagrangian (Q). It depends on the ultraviolet cutoff A 
of the effective theory in such a way as to cancel the cutoff dependence of the loop integrals 
in EQCD. 



If we use dimensional regularization in 3 — 2e spacial dimensions to cut off both ultraviolet 
and infrared divergences, then /qcd is given by the ordinary diagrammatic expansion for the 
free energy of thermal QCD. It can be calculated to next-to-next-to-leading order in g'^ by 
calculating 1-loop, 2-loop, and 3-loop diagrams without any resummation. These diagrams 
can be obtained from the results given in Ref. [^] by subtracting the effects of resummation. 



Taking the MS scale of dimensional regularization to be A and renormalizing the coupling 
constant at the scale 47rT, we obtain 

/qcd - ^ <i - ^ + — — 



+ 



IGtt 



9 [5 IQn'- 

2^^r 12 „ , A , 116 220C'(-1) 38C'(-3) 
72 log — - - 47 s V ; , s V ; 



(4) 



e ""AttT ' 5 3 C(-l) 3 C(-3) 

where 7 is Euler's constant and ({z) is the Riemann zeta function. The pole in e arises from 
an infrared divergence. The mass parameter m^j in is obtained to leading order in g"^ by 
matching the 1-loop propagator corrections for electrostatic gluons in thermal QCD to the 
propagator correction in EQCD given by the term m\AQAQ/2 in the lagrangian (^: 

N.g^T^ 



1 + e 21ogA + 2-^'(-^) 



(5) 



47rT C(-l), 

Since logarithmic ultraviolet divergences in the effective theory show up as poles in e, we 
must keep the terms of order e in (^). 

Once the parameters in £eqcd have been determined to sufficiently high accuracy, the free 
energy density F{T) can be calculated using lattice simulations of the functional integral in 
(I). If the coupling constant is sufficiently small, one might also expect to be able to calculate 
F{T) using perturbation theory. The relevant expansion parameter is gl/rriei, which is of 
order g, so the perturbation series is an expansion in powers of g{T), rather than g'^{T). 
Unfortunately, the perturbation expansion of this effective field theory breaks down due to 
infrared divergences associated with magnetostatic gluons. It is convenient to construct a 
second effective field theory containing only the magnetostatic gauge field A'^{x), which we 
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call magnetostatic QCD (MQCD). The lagrangian for this effective field theory is 



-^MQCD — -^GljGlj + (5£mqcd, (6) 



where Gij is the magnetostatic field strength with coupling constant g'^, which differs from 
(?3 by perturbative corrections. The term (5£mqcd includes all possible local gauge-invariant 
operators of dimension 5 and higher that can be constructed out of A1. Renormalization 
theory tells us that magnetostatic correlation functions at distances 1/ {gT) or larger in 
EQCD can be reproduced in MQCD to any desired accuracy by tuning the coupling constant 
(yfg and the parameters in 5£mqcd as functions of the parameters of EQCD and the ultraviolet 



cutoff A' of MQCD [T^. The partition function for thermal QCD can be expressed in terms 
of a functional integral in MQCD: 

r.(A') 



>QCD 



^T) = e-/QCD(A)T3y ^-/.qcdCA-A') (.t)3 y ^vA^^xp [- j cI^xCmqcd) ■ (7) 



In the second prefactor, /eqcdIs'^)^ can be interpreted as the coefficient of the unit operator 
in the effective lagrangian for MQCD, which was omitted in (^. Its dependence on the 
ultraviolet cutoff A' is cancelled by the A'-dependence of the loop integrals in MQCD. 

If we use dimensional regularization to cut off both the ultraviolet and the infrared diver- 
gences in the perturbation expansion for EQCD, then — /eqcd(5'^)^ is just the logarithm of 
the partition function for EQCD. It can be obtained to next-to-leading order in g by calculat- 



ing 1-loop and 2- loop diagrams in EQCD. Taking the MS scale of dimensional regularization 
to be A, the result is 



/eqcd {gTy 



47r ■' ^' I 3 ■ levrmei 



1 , A ■ 

- + 41og- + 3 

e ZrUei 



(8) 



Note that the pole in e in (H) cancels against the pole in /qcbT^, where /qcd is given in 
(^). Furthermore, after making the substitution for rriei in (1), the logarithms of A also 
cancel. The result is a logarithm of T/rriei, whose coefficient was first obtained by Toimela 



14 1 . In EQCD, the logarithm can be associated with the renormalization of the operator 



^0^0 5 which mixes with the unit operator. Thus the coefficient /eqcd of the unit operator 



in the effective lagrangian for EQCD must have a dependence on A that cancels that of the 
renormahzed operator Aq^q- 

In MQCD, perturbation expansions in the gauge couphng constant g'^ are hopelessly 
plagued with infrared divergences. Thus the functional integral in (0) can only be calculated 
by nonperturbative methods, such as lattice simulations of MQCD. Surprisingly, /mqcd can 
still be expressed as a power series in g{T). The most important momentum scale in MQCD 
is 5(3, which is of order g^T. The functional integral can therefore be expressed in the form 

where /mqcd is dimensionless. This function /mqcd has a power series expansion in powers 
of g{T). The leading term of order g^ comes from dropping the correction terms SCmqcb in 
the effective lagrangian. The resulting functional integral is 

Since a pure nonabelian gauge theory in 3 space dimensions is an ultraviolet finite theory, 
the only momentum scale is the gauge coupling constant g'^. Thus, by dimensional analysis, 
the coefficient /mqcd is a pure number that can be calculated using lattice simulations of 
the functional integral ([T(]|). Higher order corrections to /mqcd can be calculated by treating 
the terms in SCmqcb as perturbations. The leading corrections comes from the dimension- 
5 operators g'-if"'^''GijG'jk^ki and {DiGjkY{DiGjkY, which have short-distance coefficients 
proportional to g1/m\. The expectation values of these operators are of order {g'-^^^-, but 
they are ultraviolet divergent and depend on the ultraviolet cutoff A' of MQCD. This cutoff- 
dependence is cancelled by the A'-dependence of /eqcd and the parameters of >Cmqcd- Thus 
the dimension-5 operators in 5£mqcd give well-defined corrections to /mqcd that are of order 

Combining (|^) and (P), we obtain our final result for the free energy density of a non- 
abelian gauge theory: 

ES^^ = /qcd(A) + /eqcd(A,A') {gTf + /mqcd(A') (/T)^ (11) 



where g is the running couphng constant g{T). We have indicated the dependence of the 
functions on the arbitrary factorization scales A and A', which separate the momentum scales 
T, gT, and g^T. The coefficient /qcd can be calculated using ordinary perturbation theory 
in thermal QCD in the form of a power series in g'^{T), with coefficients that are polynomials 
in log(r/A). The coefficient /eqcd can be calculated using perturbation theory in EQCD, 
which gives a power series in g{T) with coefficients that are polynomial in log(A/5fT) and 
\og{gT/A'). The coefficient /mqcd can be expressed as a power series in g{T), with coeffi- 
cients that can be calculated using lattice simulations of MQCD. The functions /qcd and 
/eqcd satisfy renormalization group equations that can be used to sum up leading logarithms 
of T/A and gT/A', respectively. By choosing A of order gT and A' of order g'^T, we can 
eliminate potentially large logarithms of 1/g from the perturbative expansions of /eqcd and 
/mqcd- 

Adding (H) and (H), with g-^ = gVT and mgi given by @, we recover the result for the 
free energy of pure-glue QCD to order g^ that was recently obtained by Arnold and Zhai [0. 
The order g^ correction requires the calculation of mli to next-to-leading order in g"^ and the 
calculation of the 3-loop diagrams for the partition function of EQCD. A complete calculation 
to order g^ requires calculating the 3-loop contribution to /qcd, the 4- loop contribution to 
/eqcd, and the number /mqcd defined by (|TU|). This last number can only be obtained with 
a nonperturbative lattice calculation of the functional integral. 

The problem of the perturbative infrared catastrophe of nonabelian gauge theories is 
really one of providing a bridge between perturbative calculations of the effects of the large 
momentum scale T and the nonperturbative calculations that are required at the small 
momentum scale g^T. In this Letter, we have provided such a bridge by constructing a 
sequence of two effective field theories by first integrating out the momentum scale T and 
then integrating out the momentum scale gT. This method allows thermodynamic quantities 
like the free energy to be calculated with arbitrarily high accuracy in the high temperature 
limit. 
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